Intrinsic Spin Hall Effect in the presence of Extrinsic Spin-Orbit Scattering 
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Intrinsic and extrinsic spin Hall effects are considered together on an equal theoretical footing 
for the Rashba spin-orbit coupling in two-dimensional (2D) electron and hole systems, using the 
diagrammatic method for calculating the spin Hall conductivity. Our analytic theory for the 2D 
holes shows the expected lowest-order additive result for the spin Hall conductivity. But, the 2D 
electrons manifest a very surprising result, exhibiting a non-analyticity in the Rashba coupling 
strength a where the strictly extrinsic spin Hall conductivity (for a = 0) cannot be recovered from 
the a — > limit of the combined theory. The theoretical results are discussed in the context of 
existing experimental results. 
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The great deal of current interest in the phenomenon 
called the spin Hall effect (SHE) arises from the pos- 
sibility of controlling spin dynamics in semiconductors 
using only external electric fields. In particular, the pre- 
dicted existence of a bulk spin current transverse to the 
direction of an applied electric field (and hence the charge 
current) in doped semiconductor structures is both in- 
triguing and interesting. Recent observations of appar- 
ent spin accumulation near the edges of doped GaAs 2D 
and 3D electron Q and 2D hole jj] systems in the pres- 
ence of an applied electric field have further fuelled this 
interest as one of the possible explanations for the spin 
accumulation is that it arises from the SHE-induced bulk 
spin current although other explanations associated with 
spin precession effects at the boundaries also exist 

Spin Hall effect is traditionally theoretically discussed 
in terms of two completely distinct physical mechanisms: 
intrinsic spin Hall effect (ISHE) and extrinsic spin Hall 
effect (ESHE). ESHE, which is a rather old theoretical 
prediction, arises from the spin-orbit scattering of semi- 
conductor carriers by impurities (which is known to lead 
to spatial separation of spin-up and spin-down carriers), 
and is a solid-state analog of the well-known atomic Mott 
scattering. ESHE has recently been invoked 0, to 
explain the spin accumulation experiments in the GaAs 
electron systems although there are still some quantita- 
tive discrepancies. The intrinsic effect, ISHE, which is 
primarily responsible for the current theoretical excite- 
ment in the spin Hall effect, is an intrinsic band struc- 
ture (i.e. periodic lattice) effect, arising entirely from the 
intrinsic spin-orbit coupling in the band structure of the 
host semiconductor (i.e. GaAs) material. 

Intrinsic and extrinsic spin Hall effects have so far been 
theoretically discussed completely separately as totally 
distinct phenomena because of their completely different 
physical origins: ISHE arising from the intrinsic spin- 
orbit coupling in the semiconductor band structure and 
ESHE arising from the extrinsic impurity-induced spin- 
orbit scattering. This theoretical dichotomy is, however, 
quite odd since the two mechanisms presumably lead to 



the same observable effects, namely, a theoretical bulk 
spin Hall current, and an expreimental boundary spin ac- 
cumulation. In this Letter, we theoretically study ISHE 
and ESHE toegther on an equal footing within a single 
unified SHE formalism. Our work can therefore be con- 
sidered to be either 'intrinsic spin Hall effect in the pres- 
ence of extrinsic spin-orbit scattering' (as the title of our 
paper suggests) , or equivalently, 'extrinsic spin Hall effect 
in the presence of intrinsic Rashba spin-orbit coupling' 
(which the title of our paper could easily have been) . We 
consider the well-known and extensively studied Rashba 
model for the intrinsic spin-orbit coupling in the semicon- 
ductor carriers, restricting ourselves to 2D GaAs electron 
[2j and hole |3| systems where the Rashba coupling is ex- 
pected to be the main band structure spin-orbit coupling 
mechanism arising from the spatially imposed structural 
inversion asymmetry in the GaAs system. Our work is 
apparently the only work in the literature to theoretically 
treat both ISHE and ESHE on an equal footing within 
a unified theoretical formalism although there have been 
many studies of ISHE 0, and ESHE H Q separately. 

Our theory treating ISHE and ESHE together uses a 
minimal model with a parabolic carrier band dispersion, 
characterised by an effective mass m, with four other in- 
dependent parameters, the Rashba spin-orbit coupling 
strength a (or equivalently the Rashba spin splitting 
A), the effective extrinsic spin-orbit impurity scattering 
strength Ao, the carrier Fermi energy sf (or equivalently 
the carrier density n), and a charge transport relaxation 
time r (related to the ordinary Drude charge conductiv- 
ity through a = ne 2 T/m) completely defining the one- 
electron transport problem. We take into account both 
the side-jump (SJ) and the skew scattering (SS) contri- 
butions to the spin Hall conductivity precisely in a well- 
defined diagrammatic expansion. Below we describe our 
theory and results for 2D electrons and holes separately. 

2D electron gas with Rashba interaction. — We consider 
the single-particle Hamiltonian in the presence of both 
Rashba SO coupling and SO scattering due to impurities, 
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the disorder- free part of which is given by: 



Ho = - h ap (sm4>a x - cos^erA (1) 

2m 

where <fi = tan -1 (p y /p x ), and the disorder-dependent 
part which includes the SO scattering is given by: 



vertex correction (Fig. 1A, B, E, F), giving: 
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As shown in 0, the effect of SO extrinsic scattering can 
be taken into account by calculating the diagrams with 
corrections of the spin and charge current vertices (side- 
jump) and the SO scattering amplitude (skew scatter- 
ing). The second-quantized Hamiltonian reads 
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The spin current and charge current read 
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We proceed to calculate the ISHE conductivity a^ x using 
the Kubo-Greenwood formula, in the presence of extrin- 
sic spin Hall effect. The retarded and advanced Green 
functions are given by 
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where = p 2 /2m — ep. In the dilute impurity limit 
£pT 1, diagrams with intersection of the impurity lines 
can be neglected. We also consider the energy splitting 
between the two chiral branches A/2 = 2apF <C £f and 
neglect corrections of the order of 0(A/ef)- (These are 
all standard approximations in the context of SHE.) 

The diagrams in Fig. 1A, B, E, F and Fig. 21, J are, re- 
spectively, the diagrams for the side-jump and skew scat- 
tering contributions. Note that the Green function lines 
include the Rashba term in the Hamiltonian, Eq. 
i.e. the ESHE is now modified by ISHE. First we calcu- 
late the diagrams for the side-jump without diffuson pole 
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Diagrams in Fig. 1A, B and Fig. IE, F correspond to 
the contributions from the vertex renormalizations of the 
spin current and charge current, respectively, due to the 
anomalous SO current vertex. In the following we also 
take into account of vertex corrections due to diffuson 
poles. We have to consider the vertex corrections to two 
types of vertices, one for the charge current vertex on the 
right and the other for the spin current vertex on the left. 
For the charge current vertex, it is well known that the 
diffuson pole vertex correction leads to an exact cancel- 
lation of the spin-dependent term of the charge current 
—aa y for the Rashba model (this is the now well-accepted 
precise vanishing of the pure ISHE in the 2D Rashba 
model). Therefore, with diffuson pole vertex correction, 
the second term on the right side of Eq. J7J, which cor- 
responds to a current-current correlation of the anoma- 
lous SO current and the Rashba spin-dependent current 
— aa y , is exactly cancelled by an additional term opposite 
in sign, thus: 
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For the spin current vertex, the vertex correction r s j, 
satisfies the Bethe-Salpeter equation 
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Here N is the 2D density of states. Evaluating the first 
term under the integral on the right-hand side yields a 
term proportional to a y , which suggests that the ansatz 
r sy = jcr y , with 7 just a constant to be determined, is 
a solution of Eq. (|10|) . Straightforward calculation gives 
7 = «i?/2Ar, where vf is the Fermi velocity. Using this 
spin current vertex correction VF<J y /2Ar (Fig. 1G and 
H), we find the vertex correction to diagrams Fig. IE 
and F to be 
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FIG. 1: Diagrams for the side-jump contribution, where the 
left vertex is for y-component of the spin current and the 
right vertex is for the x-component of the charge current. The 
vertex correction due to the anomalous SO vertex is denoted 
by a dashed line connected to that vertex; while the vertex 
correction due to diffuson pole is by ladder impurity lines. 
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FIG. 2: Diagrams for the skew-scattering contribution, the 
cross denotes a correction to the Green's function Eq. @ due 
to SO scattering. 
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Finally, we have the diagrams in Fig. 20 and P which 
take into account the diffuson pole vertex corrections to 
both vertices, however, they are also the charge current 
vertex correction to the diagrams in Fig. 2M and N, 
which therefore cancels the second term in Eq. Ijl3|l , and 
< J yl +N+0+P is again given only the first term in Eq. I|13(l . 

Now we proceed to evaluate the integrals in Eqs. JSJ, 
©, Ijlll) for side-jump and the first terms on the right- 
hand side of Eqs. lfT2"|> - ljl"3l corresponding to skew scat- 
tering. Eqs. © and lfTT|) give 
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while Eqs. give 
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Now we proceed to calculate the diagrams for the skew 
scattering. Without diffuson pole vertex correction, the 
diagrams Fig. 21 and J give 
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The first sum in Eq. (|12l) corresponds to the current- 
current correlation of the spin current hk y /2m with the 
usual charge current hk x /m while the second sum corre- 
sponds to that with the Rashba spin-dependent current 
—aery. Similar to Eq. @, taking account of the vertex 
correction to the right-side charge current cancels the 
second term, so that (Jyt J+K+L is given by only the first 
term in Eq. (|T2"|) . 

The diagrams for the vertex correction to the spin cur- 



We note that these terms correspond to a combined ac- 
tion of the intrinsic and extrinsic spin Hall effect, as man- 
ifested in the appearance of both A and n^. However, 
somewhat unexpectedly, the spin current vertex correc- 
tions Fig. 1G and H cancel the side-jump contributions 
Fig. IE and F. More importantly, we find that the en- 
tire contribution to the skew scattering vanishes in the 
presence of vertex corrections as well. Now we evaluate 
the remaining Eq. 10. This gives nothing but the usual 
side-jump contribution in the purely extrinsic spin Hall 
effect: 
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We note here that the value is halved as compared to 
the case of pure ESHE since half of the side-jump 
contribution is cancelled by vertex correction. 

Without taking account of extrinsic spin Hall effect, 
the intrinsic spin Hall conductivity is (e 2 /87rft)(Ar) 2 /[l + 
(At) 2 ], which is exactly cancelled by the diffuson vertex 
correction 0. Now we have shown that, when the spin- 
orbit scattering from impurities is taken into account in 
the intrinsic spin Hall effect, there is no residual intrin- 
sic spin Hall term, nor is there any mixed-term corre- 
ponding to the combined action of the intrinsic and ex- 
trinsic spin Hall effects. Moreover, the contribution to 
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the intrinsic spin Hall conductivity from the skew scat- 
tering vanishes exactly, and only the side-jump contri- 
bution survives. Therefore, in the presence of extrinsic 
SO scattering, the total (intrinsic + extrinsic) spin Hall 
conductivity is dominated by only the extrinsic contribu- 
tion of side-jump: a^J = a SJE ^ anc j a SS _ wnere 

a xy E * s t ne side-jump result of the pure ESHE 6]. This 
result is striking - it implies that the a — > limit is non- 
analytic, i.e. SHE {a = 0) = ESHE ^ SHE (a 0) = 
[ISHE + ESHE] Q ^ - This non-analyticity, existing only 
within the Rashba model for 2D electrons, arises from 
the curious fact that as a — > 0, the vertex correction 
T sy , defined by Eq. (JTDJ, is non-zero (in fact it diverges), 
whereas for a = 0, T sy = by definition. This non- 
analyticity shows that the 2D electronic Rashba model is 
rather special. We also note that the A — > limit of the 
theory is analytic. 

2D hole gas with Rashba interaction. — We now con- 
sider a two-dimensional gas of heavy holes (with spin 
3/2) with Rashba interaction. In this case the disorder- 
free part of the Hamiltonian is given by 



H a = h ap 3 (sin30<7a; — cos34>a v ) . (17) 

2m 



The single-particle spin current and charge current oper- 
ators are 



Fig. 21, J. For side-jump, diagrams A and B yield 
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and the Green function in this case is just given by Eq. © 
with ap — > ap 3 and (f> — > 30. First we evaluate the bub- 
ble diagram without extrinsic spin-orbit interaction (i.e. 
purely intrinsic Rashba contribution), which is found to 
be o~ X y = (9e 2 /8wh)(AT) 2 /[l + (At) 2 ] (here for holes 
A/2 = 2ah 2 kp). The ladder correction to the charge 
current vertex is found to be zero. This agrees with ear- 
lier results in the clean limit. The ladder correction 
to the spin current vertex is given by Eq. I|1(J|I except 
with a extra factor of 3 attached to the bare spin current 
term. Evaluating the first term on the right hand side of 
Eq. I|10(l gives zero, because here we encounter terms with 
cos30, sin30 instead of cos(/>, sin</> in the electron case, and 
when integrated with J d<f> sin</> these terms yield zero by 
orthogonality. Therefore the ladder correction to the spin 
current vertex is zero as well, which is an expected re- 
sult since the diffuson ladder correction diagram can be 
regarded as renormalization of either the charge current 
vertex on the right or the spin current vertex on the left. 

Now that all the ladder vertex correction vanishes, we 
only have the remaining diagrams Fig. 1A, B, E, F and 
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Evaluating Fig. 21 and J for skew scattering gives 
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Evaluating the integrals, we find the second term in 
Eq. l(2U|> and also in Eq. (|2*2l are zero. The total side- 
jump contribution therefore gives: 
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and the total skew scattering contribution, within the 
short-range screened impurity assumption |6| , gives 
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It can be recognized that the combined effect of both the 
extrinsic and intrinsic spin Hall effect on heavy holes is to 
modify the extrinsic spin Hall results by a factor coming 
from the intrinsic effect. In particular, when A — > 0, 
Eqs. (|23|I - H24|) reduce to the results for extrinsic spin Hall 
effect for heavy holes. The total SHE for holes is therefore 
additive (i.e. ISHE + ESHE) in the leading order, which 
is very different from the singular result for 2D electrons. 

In summary, we emphasize that in the case of Rashba 
SO coupling for electrons, the spin Hall conductivity is 
singular at zero Rashba strength a = 0, where there will 
be a non-zero contribution coming from skew scattering 
6J. This can be traced back to the singular nature of 
the ladder vertex correction to the spin current: the so- 
lution of Eq. H10|) is zero for a = and yet becomes 
non-zero with a value of vfo~ v /2At even for an infinitesi- 
mal value of a. In the case for heavy holes, the spin Hall 
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conductivity is analytic with respect to the Rashba cou- 
pling strength and produces the expected perturbative 
result that in the leading order the net spin Hall effect 
is a sum of intrinsic and extrinsic SHE in the absence 
of each other. Finally, the implications of our theory 
for the experimental results are: (1) for 2D electrons, 
the ESHE calculated before [|| will be enhanced since 
the skew scattering contribution (with an opposite sign) 
vanishes, bringing theory and experiment in better agree- 
ment; (2) for 2D holes, the additive perturbative result 
means that both ISHE and ESHE contribute to experi- 
ment, with the ISHE being quantitatively larger. 
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